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1. Introduction 

The network of u.h.f. television stations has reached the 
stage where small pockets of unserved population are being 
provided with services from tow power relay stations. It is 
important that the sites, operating channels and specifi- 
cations of these stations are chosen so that their services 
are not impaired by interference from established or planned 
co-channel stations, and that the services of the tatter are 
not degraded by the new stations. Computer programs 
without a theoretical basis have been used for many years 
for co-channel interference (cci.) assessments, but it was 
recognised that detailed terrain information along the 
whole propagation path could give more accurate results, 
particularly for short and medium distance paths. 

The purpose of this report is to describe computer 
programs for calculating c.c.i. fields on paths up to 1000 km 
long for which profiles are obtained from a terrain data 
bank. This data, a second data bank containing transmitter 
information together with a list of receiving sites at which 
the interference is to be assessed, and the programs them- 
selves are held within the Univac 1108 computer system 
operated by the University Computing Company. The 
programs are written in Univac 1 108 Fortran V. 



2. Data banks 

2.1. Terrain data bank 

The terrain data bank consists of a large number of 
information 'units', each unit being a representative height 
and 'clutter' (building and tree) density for a small square. 
The squares are defined for three grid systems, all drawn up 
on transverse mercator projections; 



1. 



2. 
3. 



National grid of Great Britain extended to include the 

Channel Islands and the Cherbourg Peninsula (map 

origin 49°N,2°W) 

Irish grid (map origin 53-5°N, 8 W) 

'Continental' grid for N.W. Europe (map origin 

50°N, 0°) 



The data squares for grid systems 1 and 2 have side Va km 
while the continental data is on a 5 km grid since very 
detailed information about the terrain around the conti- 
nental transmitters is not required. The data is held in 
'blocks' of 400 units comprising squares of 10 km side for 
the UK and Eire and squares of 100 km side for the conti- 
nent. The extent of the terrain data contained in the data 
bank is shown in Fig. 1. The stylised continental coastline 
shown in Fig. 1 ts held within the programs for the dual 
purpose of determining the grid system in which any point 
lies and of locating the coast on paths stretching into the 
continent beyond the data limits. 



The height data itself was compiled by several 
organisations but the same rules were used to determine 
the representative height for each unit square. It is 
sufficient to note that when a square contains a definable 
feature such as a peak or valley the height of this feature is 
stored, otherwise the mean height in the square is stored. 
This scheme tends to emphasise the ruggedness of the 
terrain but this is offset by the method used to generate the 
profile (Section 4.2). For the UK and Eire, heights were 

TABLE 1 

Terrain Data Height Storage Scheme 



Height 
(ft) 


Digit Stored 


Decoded ht. 
(ft) 


Max, error 
(ft) 


(water) 

0-1 (land) 

2-6 

7-18 
19-31 
32-44 
45-55 
56-68 


1 
2 
3 
4 
5 
6 
7 
8 




1 

4 

12-5 
25 
37-5 


-\ 



1 
2 

. 6-5 


1969-1981 
1982-1993 


161 

162 


1975 
1987-5 


-/ 




1994-2012 
2013-2037 
2038-2062 


163 
164 
165 


2000 
2025 
2050 


"\ 














> 12 


3963-3987 


242 


397 


5 


^ 




3988-4024 
4025-4074 
4075-4124 


243 
244 
245 


4000 
4050 
4100 


■> 














> 25 


4375-4424 


251 


44( 


)0 


^ 
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Fig. 1 - The extent of the terrain data bank 
_2 - 



obtained from 1:25000 maps and 1 in. : 1 mile where the 
former were not available. Maps of scale up to 1:250000 
were used for the continental data. 

Clutter information w?is obtained for the UK only 
from 1 in. : 1 mile scale maps. Tree and building densities 
in each unit square were recorded on a twelve level scale. 
The occurrence of water was recorded with the tree density 
information. 

It became clear thai a total of 1-8 million units each 
with a height, tree and building density would be too large 
to handle economically. It was therefore decided to code 
the heights into 251 levels {8 bits) and reduce the clutter 
information to one bit - that is YES or NO - to enable 
four information units to be held in each 36 bit computer 
word. For economic reasons this clutter information is not 
used in tfie field calculation program. The height coding 
scheme, detailed in Table 1, utilises several features of the 
data Itself and takes advantage of the distribution of popu- 
lation with ground height in the UK. 

Maximum height accuracy is required around receiving 
sites. There are no villages in the UK higher than 1600 ft 
above mean sea level and therefore lower height accuracy 
can be tolerated above this figure. About 95% of the stored 
heights are at multiples of 25 ft and 2% are odd multiples 
of 1272 ft. The coding scheme reproduces these heights 
faithfully up to 4000 and 2000 ft respectively. The overall 
degradation introduced by the coding process is therefore 
quite small. 

The data blocks of 100 words (400 packed units) are 
stored in rows ordered from south to north through the 
three grid systems. Each block Is held In 4 sectors of a 
FASTRAND magnetic drum file and access to any block Is 
effected by means of a Fortran function 'BLOCK' which 
computes the sector address of the first word of the required 
data block from the easting and northing of the SW corner 
of the block. 

2.2. Transmitter data bank 

The transmitter data bank* consists of four randomly 
accessable FASTRAND data files which are normally 
accessed with standard Fortran read/write instructions. 

1. A Channel register contains the reference numbers of 
all stations on each channel. This file also contains a direc- 
tory for the UK relay and continental station data files. 

2. The UK main station file contains site, aerial and 
receiver location Information for the UK main stations. 
The information used within the present programs Is de- 
tailed in Table 2. 

3. The UK relay station file contains identical infor- 
mation to the main station file, but the number of receiver 
locations Is restricted to twenty. 



• The storage was devised by R.W. Lee, 
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4. The Continental station file is similar, but includes no 

receivers. 

In all the files provision has been made for separate 
aerial data for each channel, although this is frequently 
common to all channels. Storage space is also allocated to 
data used exclusively in other programs. 

3. The scope of the programs 

We shall now discuss the use of the programs in the 
planning context and describe several features of the out- 
put. 

The programs compute the protected field strength 
from an interfering source at a number of receiver sites 
associated with the service area of a station. This procedure 
is controlled by sub-program 'AREA'. The protected field 
strength (PFS) is that field which must be provided by the 
wanted station to overcome the interference from the 
unwanted source and \s the sum of the interfering field- 
strength (FS) and a 'protection ratio' (PR) which depends 
on the frequency offset characteristics of the two trans- 
mitters and the aerial pattern {polarisation and directivity) 
of the receiving aerial. 

Two types of receivers are recognised for c.ci. 
purposes: domestic receivers, at height 10 m (32-8 ft) 
above ground level, and receivers on re-broadcasting links 
(r.b.l.) which may be at any height. The latter are all 
located at transmitter sites and their data is accessed from 
the transmitter data bank. The two receiver types have 
different stylised aerial directivity patterns for protection 
ratio calculation. For all receivers the median field strength 
iFS,.) at the aerial is computed together with the field 
strength exceeded for 5% [FS^] of the time for domestic 
sites and that exceeded for 1% (FS^) of the time for r.b.l. 
sites. The PFS is calculated for the worst case; thus for a 
domestic site 

PFS = MAX (KSjo + 10, F5j )+PR dB /iV/m (1^) 

and for an r.b.l. 

PFS=MAX{FS^^ + W,FS^)+PR dBMV/m db) 

The 10 dB factor added to FS^^^ is an allowance for 
persistent interference. 

A number of operating versions are available,, but 
the two more important are: 

1. Specify one station and compute the c.ci. from or to 
all stations operating on the channels of the specified 
station. The two types of prediction are 'forward', 
i.e. the c.ci. at the receiver locations of the specified 
station from all co-channel stations, and 'reverse', i.e. 
the c.ci. from the specified station to the receiver 
locations of the co-channel stations. 

2. Specify two stations and compute the c.ci. from the 
first to the second on a single channel. More 
receiver locations can be used in this version and a 
graphical output may be provided. 
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Fig. 2- Ouput from version 2 prediction 
[h) plotted output 
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FILES 
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Fig. 3 - Simplified flowchart for version 1 control program 
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TABLE 2 
Transmitter Data 



Item 



transmitter number 

transmitter name 

responsibility 

active deflector indicator 

program status 

tx n.g.r. 

latitude 

longitude 

site height (ft a.m.s.l.) 

polarisation 

date of data 

number of r.b.i.'s 

number of domestic receivers 

r.b-i. numbers 

domestic rx. n.g.r.'s 
FOR EACH CHANNEL: 

channel number 

offset [^h line) 

precision indication 

tx. aerial height (ft a.g.l.) 

rx. aerial height (ft a.g.l.) 

program source 

max. e.r.p. (kW) 

aperture (wavelengths) 

beam tilt {° below horizontal) at 20 

intervals 

h.r.p. (dB below max.) at 10° intervals 



Type 



A 
A 
A 
A 



No. of Words 



1 

3 

1 
1 

1 
1 
1 

1 
1 
1 
1 
1 each 

1 each 

1 
1 



Notes 



2 letter, 6 figure n.g.r. packed into one 
word 



Total number of receivers: 
60 for main stations 
20 for relay stations 
for continental stations 
packed as for tx. n.g.r. 



condition allowed: precision, 
stabilised, unstabilised 

> for re-broadcast transmitters 



tilt for three directions packed into 
each word 

coded figure for six directions packed 
into each word 



I : integer, F : floating point, A : character format 



A user may operate the versions in several ways. For 
example he might run version 1 'forward' on a single 
channel and find one particular station gives a higher field 
in the service area of the specified station than he had 
anticipated. He might then run version 2 with the parti- 
cular station causing the interference and include more 
receiver locations. The user might, at this stage, change the 
parameters of the specified station and re-run. Finally to 
obtain a complete assessment of the station he would run 
version 1 on all channels of the specified station in both 
'forward' and 'reverse' directions. 

Because of its simplicity we shall describe the output 
of version 2 first. As the example in Fig. 2 shows, this con- 
sists of details of the two transmitters, followed by the 
receiver details: name (for r.b.i.'s), grid reference, FS^^, 
PFS and distances and bearings to the transmitters. The 
'W beside the PFS value indicates thati'FS arises from the 
median field value, rather than the short percentage time 
value. The accompanying plot (Fig. 2(b)) gives the values 
of PFS for the domestic receivers on 1:25000 or 1 in. : 1 
mile scale. 



A simplified flow chart of the version 1 control 
program is shown in Fig. 3. The output for this version is 
somewhat more complicated than that of version 2, con- 
sisting of details of the specified station and all co-channel 
stations.and the details of the field at the receiver locations, 
as shown in Figs. 4 ('forward') and Fig. 5 {'reversals'). The 
first part of the output (Fig. 4{fl) and Fig. 5(a)) consists of 
the specified station name, channels and offsets, status of 
each channel and date of last data change. This is followed 
by similar details of the co-channel sources. The column 
PR with each channel gives the protection ratio (dB) on 
that channel for persistent interference, taking into account 
offset conditions at the two transmitters. A non shared 



channel has 



in ii\\ePR column. 



In order to eliminate a large number of unnecessary 
calculations (from low power distant sources, for example) 
each transmitter is tested with a procedure to determine 
whether or not its influence on domestic receivers is likely 
to be significant. 

The 'exclusion procedure' is based on the path 
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between the unwanted and wanted transmitters. There are 
four tests that the path has to pass before the unwanted 
transmitter is judged to be of importance. The number of 
the test at which the exclusion occurs is printed in column 
EXC/PROC on the listed outputs. The four tests are 

1{a} path distance >1000 km (no distance or value for 
PFS ]s printed) 

Kb) protected field strength at the wanted transmitter site 
<x dB /iV/m, assuming the path to be a sea path. 
This calculation of protected field strength is based on 
the field strength exceeded for 5% of the time and the 
protection ratio does not include receiver aerial direc- 
tivity. 

2(a) wanted transmitter in N. Ireland, unwanted transmitter 
in Continental Europe. The value of PFS listed for a 
a path excluded on this basis is that from 1 (b). 

2(b) protected field strength at the wanted transmitter site 
<jc dB AtV/m, for the path extracted from the terrain 
data bank. The field strength calculation is for 5% 
time, and the protection ratio does not include aerial 
directivity. 

The basic value of x is 40 dB nV/m when the wanted trans- 
mitter is a high power station with a large sen/ice area and 
50dB/jV/m for a low power relay station. A further 10 dB 
is added to this value if the distance between the two trans- 
mitters is>200 km. 

The exclusion procedure Is normally successful in 
removing about 80% of the interfering sources. 

The listing of the conditions at the receiver locations 
for 'reversals' (i.e. influence of the specified station) shown 
in Fig. 5(6) is similar to that already described for version 2 
output wilh the addition of the number of the channel 
which produces the highest protected field. The frequency 
used in the field strength calculation is appropriate to the 
arithmetic mean of the common channels. 

The listing for 'forward' prediction {i.e. the influence 
on the specified station} is shown in Fig. Mb) and 4(c) for 
domestic and r.b.l. sites respectively. The components of 
the output which are given for those stations not excluded 
are compatible with those described above. For r.b.l. 
sites, the ou^ut is not compressed across the page because 
each transmitter is handled individually for each r.b.l. for 
exclusion purposes. The exclusion is at 70 dB iiV/m for 
PFS at the r.b.l. site, including aerial protection, for the 
path assumed to be sea. The calculation is for 5% time. 
The printed list includes only those transmitters not 
excluded on this basis or on conditions 1(a) or 2(a) as 
above. (The r.b.l. site replaces the wanted transmitter in 
the above nomenclature.) 



4. The control sub-program (AREA) 

The purpose of this sub-program is to control the calcu- 
lation of protected field strength from an unwanted trans- 



mitter at the receiver locations of a wanted transmitter. A 
simplified flow-chart of its operation is shown in Fig. 6. 
Much of the sub-program is involved in sorting the un- 
wanted transmitter to receiver paths into categories and 
controlling terrain data handling and generating path pro- 
files as a series of heights and distances (subroutines 
TRANAR, SPATH and PATH). The final part of the sub- 
program controls the field strength calculation (GRIDFS), 
radiation pattern correction (AECORR), protection ratio 
calculation (PROTCT) and stores the results prior to print- 
ing. The named subroutines will be described later. 

It was found to be uneconomic to access each path 
from the data-bank directly so terrain data is always trans- 
ferred into core-store in bulk (TRANAR) before the paths 
are extracted. However the area of core available for this 
data is limited to about 20,000 words, so that long and 
short paths are handled in slightly different ways. For 
unwanted transmitters within 150 km of a group of receivers 
the data for all the paths can generally be transferred and 
held in core whilst the profiles are extracted (PATH) and 
protected fields calculated for each path in turn. For 
longer paths, controlled by SPATH, only sufficient data for 
150 km lengths is held atone time, so each path is built up 
in 150 km sections (TRANAR and PATH) and the calcu- 
lations are completed before proceeding to the next path. 

4.1. Transfer of terrain data to core-store (TRANAR) 

The procedure is illustrated in Fig. 7. First, the 
extremities of the area enclosing the receivers are deter- 
mined and a rectangular box is constructed to, enclose the 
data blocks (10 x 10 km squares for grid systems .1 and 2, 
100 X 100 km square for grid system 3) containing the 
receiver area. The transmitter point is surrounded by a box 
of side 1 km (10 km for grid system 3) to allow for data 
interpolation along the path, and the data blocks enclosed 
are also determined. The data blocks to be transferred to 
core are determined by the eastings of the intersections of 
lines AA' and BB' with the data rows crossed. A simple 
directory is then established. This contains the most 
westerly (IW) and most southerly (IS) grid references of the 
required data blocks, and for each data row i, a count of the 
first square in that row lC(i) and its westing lD(i) relative 
to IW. IC(i) for an extra row beyond the data is also 
required in order to record the number of squares in the 
last row. The director/- for the data of Fig. 7 shown in 
Tables, performs the same function for the transferred data 
in core as does the function BLOCK for the terrain data 
bank. 

TABLE 3 
Director/ for the data transferred from drum to core (Fig. 7) 
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Fig. 7 - To illustrals terrain data transfer to core store 

Having determined the data blocks required, the 
program transfers the data row by row. If a data block is 
not in the data bank, a dummy block containing 'sea' is 
inserted into the core data. 

A two stage transfer requiring two distinguishable 
areas of core is necessary when paths cross from one grid 
system to another. 

4.2. Path profile generation (PATH) 

A path profile can be interpolated from terrain data 
in several ways. Methods examined within the present 
context were 

(a) ortfiogonal projection of the data point on to the path 




grid tieights assumed to be 
at centre of grid squares 
projected orlhogonoMy onto 
profile line to determine 
distances 



0-5 km 



In method (a), illustrated in Fig. 8, the height stored for 
each data square through which the path passes is projected 
orthogonally from the centre of the square on to the path. 
The resulting profile therefore retains the exaggerated 
ruggedness of the stored data. The profile points have non 
uniform separations. 

Method lb}, on the other hand forms a profile with 
uniform or specifiable separation between the points, and, 
as illustrated in Fig. 9, linearly interpolates a height from 
the stored heights for the four data squares surrounding the 
point. The height at P is given by 

- {(h^ -h^}x + h^]y (2} 

A surface defined in this way has the appearance of a 
'twisted plane' and any section across it parallel to eitfier 
grid axis is a straight tine. As this height is the result of 
interpolation the profile is somewhat less rugged than that 
obtained from method (a). Method (c) uses row, column 
or diagonal interpolation along the profile, according to 
rules which depend on the angle of the path line with the 
grid. The method is complicated and can result in different 
heights for the same point (depending on the interpolation 
direction) and does not produce a uniform distribution of 
profile points. 

Method (b) was chosen because it 'softens' the data, 
produces a unique profile and is economic in computing 
time and program storage. 

The profile heights are generated at intervals depend- 
ing on the distance from the terminals. 



Thus, 



Fig. 8 - Profile generation by ortiiogonal projection 
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'^i+l =^i + ^Anc 



where 



60<rfi'<120 

10<di'< 60 

0<d{< 10 

0<d< 2 



^inc=4-0 
^inc = 2-0 



:th 



where d-^ is the distance of the i point from the receiver 

(km) 
rfj+l is the distance of the (i+D* point from the 

receiver (km) 
d-l is the distance between the i* point and the 

nearer terminal (km) 
rfjric is the distance interval (km) 

No quantitative comparison between true (map) and data 
bank profiles has been attempted, but the accuracy is 
obviously reflected in the comparison between measured 
and predicted fields described in the Appendix. 

Profiles of paths <150 km and each 150 km section 
of longer paths are generated as arrays of heights and 
distances in subroutine PATH, together with details of 
coastal intersections. A stretch of sea is recognised as such 
if it occupies more than 5 km along the path. 

4.3. Paths longer than 150 km or extending beyond the 
terrain data bank (SPATH) 

Profiles for these paths are generated in sections of 
length 150 km because of limitations of terrain data core- 
storage. The control splits each patfi into sections and, 
using TRANAR and PATH, extracts the profile for each 
section in turn. 

Paths with the transmitter beyond the limits of the 
data bank have to be handled in a different and less accurate 
way. These paths are never shorter than about 200 km 
and generally contain considerable stretches of sea. The 
procedure used in 'SPATH' is to find the intersection X 
between the path and the stylised continental coastline 
shown in Fig. 1. X is then used as a temporary terminal to 
generate the path profile from the receiver end using 
'TRANAR' and 'PATH' in ttie normal way. The profile 
is completed by adding points at 20 km intervals from X 
to the transmitter at the given transmitter ground height. 

4.4. Field-strength calculation (GRIDFS) 

The philosophy of the field calculation has been 
described elsewhere. Briefly, it involves computing a path 
diffraction loss ^q, as a spherical loss/l^ (for sea paths) or 
as an interpolation between a multiple knife-edge loss Ay 
and a loss calculated over a stylised surface (wedge, A^ or 
composite cylinder, A^) for land paths. Except on sea 
paths for field levels exceeded for 1% and 5% of the time, a 
loss -4 -y- is computed from the path angular distance. This 
is to allow for fields resulting from tropospheric irregu- 
larities. The basic path loss is the minimum of /4q and /I-j- 
with the sea path exceptions above, for which loss is com- 
puted on a distance basis. The path loss is completed by 
adding a 'clutter loss', Aq_^. A simplified flow chart of 
subroutine GRIDFS is shown in Fig, 10. 



The input data to the subroutine consists of the path 
profile — heights above mean sea level, distances and 
coastal intersections and the wavelength. The first opera- 
tion is to transform the profile heights into a cartesian 
system, based on a number of effective earth radii. These 
radii are obtained from the true radius, taken to be 6367 
kilometres, by multiplying the latter by a value depending 
on the percentage time for which the calculation is being 
made, and on whether the terrain is land or sea. The 
multipliers used are those obtained in the optimisation 
process described in the Appendix. 

An imaginary string is stretched between the two 
terminals over the transformed profile to determine which 
points, called running edges, touch the string. The calcu- 
lation procedure depends on the number and distribution 
of these running edges. But first the effect of clutter is 
assessed by computing the occupancy by buildings and trees 
within 2 kilometres of the receiver in a zone surrounding 
the 'string'. This zone is smaller than the first Fresnel zone. 
The method is based on that given in Reference 4. The 
density of 0-5 was taken for the land value and when the 
path crosses over the sea. The height distribution was 
taken as that given tor trees in Reference 4. The constant 
value of 3 dB is included in the final clutter loss, A^y^, 
except for cases where the coasf is within 2 kilometres of 
the receiver, or the receiver is an r.b.l. 

The running edges are tested against three 'grouping' 
criteria to determine whether or not certain edges are taken 
as a single crest: 



d: 



i+i 



d, < 0-5 km 



(b) 



dSD-d-^,,) 
d-.^^^D-d.^) 



> 0-9095 



(0 



with D =^ total path length (km) 

and d-^, d-^+-^ are a pair of adjacent running edges (km) 

as (a) and (b) above but instead of running edges three 
points are considered: the receiver horizon, the 
junction point of transmitter and receiver horizon ray 
lines, and the transmitter horizon. 



If a set of points requires grouping on any of these bases, a 
virtual edge is formed using the first and last edges of the 
group. 

The procedure now depends on the number of 
running edges after grouping and a 'flag' which indicates if 
any of the edges are formed by the sea. The following 
categories are recognised: 

I No edges (line-of-sight) 

Here, the subroutine DBMAXK determines the path 
profile point which gives the maximum Fresnel parameter 
V (defined in Reference 3). Having found this point three 
conditions are recognised: 

(i) [-'<— 0-78 : Diffraction toss taken as zero 
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(ii) Pr>— 0-78 with land forming maximum V point : 
a wedge loss A^ is calculated. The maximum V 
point is taken as the wedge apex. The wedge surfaces 
are defined by the parts of the profile stretching from 
the terminals to distances short of the apex given by 
grouping criteria (a) and (b). The surfaces formed 
touch those points in this profile which produce the 
maximum internal wedge angle. A restriction is that 
t he we dge surface shall be kept below the aerials by 
V45cf' feet, wtiere d' kilometres is the distance from 
the particular terminal to the wedge apex. A four-ray 
approach, explained in the Appendix, is then used to 
complete the wedge loss calculation. 

(iii) K>— 0-78 with sea forming maximum V point : A 
spherical diffraction formula is used to calculate 
fields exceeded for 50% of the time. The radius of 
the sphere is taken to be 8490 km, i.e. h the 
physical radius. To obtain fields exceeded for 1% 
and 5% of the time tropospheric ducts of limited dis- 
tance are assumed. The loss is taken as: 



A^ (1%)= MAX (0, 94-2 Log^^D - 2441 
A^ (5%}= MAX (0, 82 Log,^/)- 179) 



(3) 



This formulation allows free-space propagation for 
383 and 153 km respectively. 

II One edge 

Here, two major conditions are recognised: 

(i) edge on land : a wedge calculation is made as in I {ii}, 
except that the wedge surfaces are defined by the 
parts of the profile short of the horizons by the 
stated amount. 

(ii) edge formed by the sea : the loss, .4^ depends on the 
time percentage for which the field is required. To 
obtain field levels exceeded for 50% of the time, the 
lossiscalculated from aspherical diffraction formula. 
The levels exceeded for 5% and 1% of the time are 
determined via Equations (3). 

I i I Two or more edges 

A^^ h calculated for a series of knife-edges by the 
method of Reference 3. The profile is then stylised as a 
surface of four right-circular cylinders forming a single 
convex surface in which there is no discontinuity of slope. 
The radii of the inner cylinders are determined by the 
position of the horizon and the direction of the horizon 
rays, and the outer cylinders depend on a similar part of the 
profile to that selected for the wedge stylisation. The 
cylinders are chosen to pass through profile points which 
produce the maximum radii, subject to the same restriction 
that was applied to the wedge surface. The composite 
cylindrical loss, A^, is then derived from a method due to 
Vogler, described in Reference 3. The diffraction loss is 
computed from: 



if no sea edges flagged, and 

A^ = 0-24A^^ + 0-6 MIN (A^. 42 Log,^ D] (5) 

if sea edges are flagged. 

For both categories 2 and 3 above, but with the 
exception of sea path fields exceeded for 5% and 1% of the 
time, a determination is made of the angle 6 between the 
terminal horizon rays, if this angle exceeds 2 mr a scatter 
calculation is made using: 

/l-r-MAX{10+37Log,gf)^„28+ 10 Log, „ d^^) [6] 

A justification of this form is given in the Appendix. 

The final pa^ attenuation relative to free space is 
given by: 

A - MIN {Aj.Aq) + A^^_ (7) 

and the field for 1 kW e.r.p. is given by: 

/•■= 106'9- 20Log,^ D A dBA(V/m (8) 

4.5. Radiation pattern correction (AECORR) 

Account of the power radiation pattern of the trans- 
mitting aerial is taken by linear interpolation in the 
azimuthal direction between the values of relative radiated 
power (dB) stored at 10° intervals, and calculating the 
vertical pattern from a modified sine function, based on 
beam tilt and aerial aperture. Finally 



FS = F-i-P- H - V dBMV/m 



(9) 



A^ = 0-44/1 1, + 0-38 



UN U^, 42Log,„/)l 



(4) 



where/* is the maximum e.r.p. in dB relative to 1 kW, 
H is the loss in dB in the azimuthal plane 
V is the loss in dB in the vertical direction 

4.6. Protection ratio calculation (PROTCT) 

The two components of the protection ratio have 
already been mentioned. The first component depends on 
the frequency offset of the two transmitters and the per- 
centage time for which the calculation is being made. 
10 dB extra protection is required to overcome persistant 
interference, compared with that required to overcome that 
present for a short percentage of the time, quoted in 
Reference 1. 

The receiving aerial protection depends upon the 
type of receiver (domestic or r.b.l.) being considered. The 
patterns vary in azimuth only, and are used with the 
assumption that the receiving aerial is always directed 
towards the wanted transmitter. 

Equation (1) is used to convert field-strength to 
protected field-strength. 
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5. Input data 

The programs are designed to operate with the mini- 
tnutn of input data at the time of the run. Thus for 
version 1, for example, it is necessary only to specify the 
number of the station under consideration, the direction in 
which the calculation is to be made ('forward' or 'reverse'), 
and a code to indicate whether the transmitter data is to be 
accessed from the transmitter data bank or Input on furtiier 
cards. A second code indicates whether or not the receiver 
data held with the transmitter information is to be supple- 
mented by further card data. Similarly for version 2, both 
stations may be accessed from the data bank or input from 
cards, and additional information consists of the channel 
under consideration and a map scale code if the results are 
to be presented graphically. 

. The format of the cards containing transmitter Infor- 
mation is identical to that used In setting up and updating 
the stored transmitter data bank. A set of cards for each 
station Is available. Thus, the effect of a small change, such 
as the offset condition,^ of a stored station can be quickly 
assessed by running the program with the data for the 
station input from a set of cards Including the modification. 



6. Conclusions 

The programs described In the report enable cc.l. 
levels to be predicted in the UK from stations in the British 
Isles and N.W. Europe with acceptable accuracy and reason- 
able economy. The use of randomly accessable terrain and 
transmitter data banks held on a mass storage drum enable 
predictions to be made very rapidly. 

The system Is deficient In a number of ways. The 
true variation of field within 14 square km can be quite 
large and the figure quoted may not be truly representative 
because: 

(i) the height stored for the square may not be represen- 
tative of the heights of the households in the square 

(ii) the profile within a short distance of the receiver may 
not be a faithful representation of the true (map) 
profile, because of coarseness of the stored data. 

To a certain extent these deficiencies can be overcome by 
examining the relative fields in adjacent squares. 

Increased accuracy could undoubtedly by obtained 
with extremely detailed profile data. This lack of detail 
made it impossible to calculate losses based upon the 
curvature of individual hill tops."^ The requirement for 
this would be contour heights at 5 m intervals mapped 
with a horizontal accuracy of better than ±10 m. This 
requirement Is beyond the capability of current maps. 
Accuracy could also be increased with more clutter data but 
only at the cost of a considerable increase in data storage 
and handling charges. This deficiency becomes important 



when predicting the signals from low power fill-in stations 
with a small coverage. Thus in these cases calculations are 
made with the programs described in Reference 4. 

The calculation system could be improved if more 
were known about propagation in the region where the pure 
diffraction and tropospheric mechanism give roughly equal 
fields. We do not know to what extent terrain on either 
side of the profile line can change the field In any given 
situation. The tropospheric mechanisms used In the pro- 
gram are rather cmde, and probably greater accuracy could 
be achieved by allowing more phenomena, e.g. trapped 
propagation over low tying land, to contribute to the field. 

Nevertheless the field-strength calculation has the 
philosophical advantage over earlier methods that the cal- 
culation method used on any particular path is based on the 
dominant propagation mechanism on that path, as far as it 
can be known. The distribution of differences between 
measurement and prediction indicates that the calculation 
system is more or less uniformly successful over the wide 
range of path attenuations encountered in practice. 
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Appendix 
Minimisation of Differences Between Calculation and Measurement 



The system used to calculate attenuation relative to 
free space has a theoretical basis, but empiricism is also 
necessary to overcome imperfections in applying the 
theories. The Appendix outlines the measurement projects 
used and the processes adopted to develop a calculation 
system with a small difference from these measu.rements. 

The field prediction program is required to give 
variations in both time and space, but not for precise times 
or precise points in space. The time variations required are 
field levels exceeded for 50, 5 and 1% of the broadcast 
time, and the primary information required about location 
is the mean value of the fields at households with in areas of 
% sq. km. 

At short distances, where time variability is small, 
each field value used was the mean of many measurements 
in the relevant area. Most of the field measurements were 
made in Guildford and the Potteries. In Guildford the field 
strength was measured from transmitters at Reigate, Oxford, 
Guildford (relay). Crystal Palace and a small transmitter 
with low aerial specially installed on the Crystal Palace mast 
to compare with the main transmitter. In the Potteries the 
signals were received from Winter Hill, Sutton Coldfield, 
Brierley Hill and Fenton. Values also came from measure- 
ments taken from the Rowridge transmitter just beyond its 
service area. The total number of areas of this type was 
approximately 500. 

Evidence on time variability came from experiments 
conducted by several organisations. Only those paths 
which were measured for sufficiently long time periods 
were accepted as evidence. The measurements at these 
fixed sites were not always typical of the area, so most 
experimenters obtained site variation factors (S.V.F.) by 
taking a series of measurements in a district and relating 
tiiem to a fixed site. It was decided not to use these factors 
because the areas they represented were poorly defined and 
were not specific to the present requirement. 

If it is assumed the fixed site does not suffer from 
muitipath effects and the ground profile changes are not 
too abrupt, then the only features liable to produce fields 
untypical of the surrounding % sq. km are the buildings and 
trees. For this reason the clutter loss was determined 
manually using data derived by visits to the sites. Details 
of the patfis used for time variability are given in Tables 4, 
5 and 6. 

Propagation over sea is substantially different from 
propagation over land, so paths were categorised with this 
in mind. This categorisation is aided by imagining a string 
stretched across the profile from transmitting to receiving 
aerials. Thus: 

(1) String touches only land and no sea stretches in the 
path exceed 5 km in length. 



(2) String touches only land, and sea stretches in path 
exceed 5 km in length. 

(3) String touches land and sea. 

(4) String touches only sea. 

If none of the profile points is touched by the string 
but a significant loss is given by terrain below it, the cate- 
gory used is 1 or 2 if tiie largest intrusion is land and if it is 
sea category 4 is taken. If no intrusion into the zone by the 
earth takes place categorisation is irrelevant. 

In the tables the first categon/ is labelled 'land', the 
second and third are labelled 'mixed' and the fourth 'sea'. 
The categorisation of a path may change if the effective 
curvature of the earth is changes, a concept covered below. 

For all categories and time percentages the calculation 
of attenuation relative to free space {A in dB) is obtained 
by: 



A=A^^_ + M\N (Aj.Ao'. 



;io) 



where^QL " '°^^ caused by clutter (trees and buildings) 

Aj = attenuation relative to free space of field 

received due to irregularities of troposphere 

Aq = loss caused by diffraction of earth's surface. 

The values obtained depend on tfie radius assumed for 
the earth. The time variations are achieved by assigning 
varying effective radii to the earth. The radius assumed on 
both land and sea for the mean field level is four thirds the 
physical radius, but for other time percentages the effective 
radii are determined empirically. 

Prediction of Aq\_ is based on a method described in 
a companion report on service area prediction. However, 
in this case several simplifications have been made for 
economic reasons. The loss is given by: 



'4cL = ^ +'^,VCg + X3C 



where 



Cg = 



=E E 



(//l 



f(>' 



;n: 



;i2} 



The summation of distance, over i, starts at 55 metres and 
then from 200 metres to the coast or 2 kilometres (which- 
ever is shorter) in 200 metre steps, ^j is the radius of the 
volume being considered at the distance in question, with 
centre at //j. The summation of height (j) is from the 
ground or bottom of ttie volume, (whichever is higher) to 
the top of the volume in steps of 1 metre, y-^ is the height 
of the strip above ground level. The function 'f is given by 
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Fig. 1 1 - Probability, f, of clu tter exceeding a given heigh t 

Fig. 11. The values x,, x^ and x^ were determined 
empirically. 

The parameter used to calculate A-^ is the angle, 6. 
between the horizon line from the transmitter and the one 
from the receiver, the functional relationship being deter- 
mined empirically. The evidence for this work was selected 
by examining the original field recording charts for the 
typical fast fading. The selected points are shown in Fig. 
12 as median time attenuation relative to free space plotted 
against log Q^^. Correction was made for clutter losses 
before the values were plotted on the graph. From this the 
decision was made to calculate the tropospheric loss as, 

^^ = MAX (10 + 37 log 0^,, 28+ 10 log 0^^) (13} 
This gives the lines drawn in Fig. 12. 

A^ is calculated for categories 1, 2 and 3 by means of 
various stylisations of the shape of the terrain profile i.e. 
knife-edges, wedges and cylinders. 

If after the grouping process, described in Section 4.4, 



p 40- 




Fig. 12 - l^edian losses occurring due to tropospheric 

irregularities for paths categorised into three types 

O sea X mixed A land 



only one edge exists between transmitter and receiver, the 
wedge stylisation is adopted. The value of A^ is then 
given by: 

(14) 



where 



F(i',} = 



1-i 



r 77 



,2 ' df 



(15) 



and the Vj are the variables associated with the terminals 
and their images as described in Reference 3. The values of 
R^ and R^ are simulated reflection coefficients, for each 
side of the wedge, which are taken as functions of the 
average differences in height h^ and h^ (feet) between the 
profile and the assumed wedge. The relationship between 
the h and R values is taken to be 



R = 



-0-99 



1 + x^/i + x^h'^ 



:i6) 



where x_j and x were determined empirically. 

If after grouping there remains more than one edge 
between transmitter and receiver, multiple knife-edge and 
cylindrical calculations are made and the value of Aq 
determined by interpolation between the losses given by 
these two stylisations. Thus: 



A^ = x^A^^ + X, MIN(y4^, 42 iog,^^) 



;i7i 



The term in which D (the path length in km) is used was 
included because certain data anomolies produced large 
overestimates oi A^. 

In category 4 the value of A^ is calculated directly 
assuming a simple spherical profile. 
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fig. 13- Comparison between predicted and measured path attenuations for paths on which small time variation was assumed 



Values Xj to x^ for categories 1 and 2 were deter- 
mined in a computer optimisation technique, based on a 
method by Dickinson, which minimised the sum of the 
squared difference between calculation and measurement 
over all cases. Approximately 550 paths were used, with 
the result that: 



[18) 



These results gave differences between calculation and 
measurements having a standard deviation of 7 dB with a 



X, - 


3-0 


X. = 


5-93 


x^ = 


-0-475 


^'i = 


0-0045 


Xr = 


3-04 X 10"^ 


''r = 


0-44 


X, = 


0-38 



spread from 
Fig. 13. 



■18 to 19 dB. These results are plotted in 



In category 3 the path loss is expected to be nearer to 
the cylindrical value than to the knife-edge result. The 
empirical evidence for this was that on such paths Equation 
17 underestimated Aq to such an extent that diffrac- 
tion rather than tropospheric scatter became the dominant 
mechanism. Field charts for these paths showed that the 
latter mechanism was prevelant, and the anomaly was 
avoided by replacingx and x , for such paths, by: 



X, = 0-24 



x„ = 0-6 



[191 



The determination of effective radii for the 5 and 1% 
time values was carried out in a heuristic manner. The 
distinction between categories 1 and 2 is necessary because 
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TABLE 7 



Effective Eartfi f^adius Multiplier, k 



% time 


land 


sea 


50 
5 

1 


1-33 
2-2 

4-5 


1-33 
10-0 
25-0 



the land and sea sections {in 2) were allowed different 
effective radii. (The sections join with a common tangent.) 
The multiplying factors are given in Table 7. The calcu- 



84 
80 
76 
72 
68 

64 
60 
56 
52 
48 
44 
40 
36 
32 
28 
24 
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12 
8 
4 
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lation of A is then the same as for 50% time except paths in 
category 4. Here tropospheric ducts are liable to influence 
propagation, and the path loss is taken to be of the form 
given in Reference 1. For 5% time this is: 

A =/lcL + MAXJO, (82 togD- 179)[ (20) 

and for 1% time:^ 

A=A^^+ MAXJO, (94-2 log D- 244)} (21} 

In the case of long sea paths it Is not reasonable to 
suppose that the high multiplying factors, given in Table 7, 
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12 16 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80 84 
measured attenuation relative to free space, dB 



Fig. 14 - Comparison between predicted and measured path attenuations for field levels exceeded for specified time 

percentages for paths measured over long periods of time 
1% X 5% • 50% 
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are relevant. Thus a reduction is made to these values if interpolation. 

the sea sections total more than 153 km on a 5% time 

calculation, and 383 km on a 1% time calculation. It is The differences between calculation and measurement 

assumed that if the sea sections total more than 1000 km produced by these techniques are shown in the last three 

the equivalent land multiplying factor is a reasonable columns of Tables 4, 5 and 6. At the foot of each column 

estimate. Then for a length between the 153/383 and the is the mean error and the standard deviation. These 

1000 km distances the multiplying factor is determined by results are also presented in Fig. 14. 
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